Spatial-Spectral Vortex Solitons in Quadratic Lattices 
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We predict the existence of spatial-spectral vortex solitons in one-dimensional periodic waveguide arrays with 
quadratic nonlinear response. In such vortices the energy flow forms a closed loop through the simultaneous 
effects of phase gradients at the fundamental frequency and second-harmonic fields, and the parametric 
frequency conversion between the spectral components. The linear stability analysis shows that such modes are 
stable in a broad parameter region. © 2009 Optical Society of America 
OCIS codes: 190.0190, 190.4420, 190.6135 
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Manipulation of light beams and pulses in nonlinear 
photonic lattices or waveguide arrays is attracting in- 
creasing attention, due to the potential to control spatial 
beam shaping combined with manipulation of temporal 
and spectral characteristics [1]. In particular, photonic 
lattices created in a medium with quadratic nonlinear- 
ity can facilitate ultra-fast all-optical switching through 
parametric wave mixing between the fundamental and 
second-harmonic waves [2] . Various approaches to beam 
manipulation reply on the special features of localized 
modes in the form of discrete or lattice solitons [1,3-6]. 

In this Letter, we predict the appearance of a differ- 
ent type of discrete quadratic solitons in one-dimensional 
lattices. In contrast to previously known solitons, such 
localized states exhibit directional power flows between 
the lattice sites, that are compensated through para- 
metric conversion between the fundamental and second- 
harmonic waves. Since the power flow loop is closed in 
spatial-spectral domain, we refer to these solutions as 
spatial-spectral vortices. 

The spatial light dynamics in a one-dimensional pho- 
tonic lattice created in a medium with quadratic non- 
linearity can be modelled by a set of coupled equa- 
tions for the amplitudes of fundamental frequency (FF) 
and second-harmonic (SH) modes of individual waveg- 
uides [1], which can be written in the normalized form: 
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where z is the normalized propagation distance along 
the waveguides, An and _B„ arc the normalized FF and 
SH mode amplitudes in the nth waveguide, respectively. 
The coefficients ci_2 define the linear coupling between 
the guided modes at the corresponding frequency compo- 
nents. Parameter (3 defines the phase mismatch between 
the FF and SH modes. The strength of quadratic non- 
linearity is characterized by the coefficient 7, and with 
no loss of generality it can be scaled to unity (7 = 1). 
In order to study vortex states, it is important to con- 



sider the mechanisms of power flows. The power den- 
sity at the lattice site number n can be defined as 
and |-B„P at the FF and SH spectral components, re- 
spectively. Whereas it can be shown that the overall 
power is conserved, / — (^|^n|^ + |^n|^^ = const, 
the power density at individual lattice sites can change 
due to (i) mode coupling with the neighboring lattice 
sites of the same frequency component and (ii) paramet- 
ric frequency conversion. Indeed following from Eq. (fT|). 
we have d\An ^ jdz = Ji{n ~ l,n) ~ Ji{n,n + 1) — Jp(n), 
and d\Bn\'^/dz = J-iin — l,7i) — J2{n,n + 1) + Jp{n), 
where Ji(n, = 2ciIm(A* A„+i) and J2{n,n + 1) = 

2c2lni(i?*i?„+i) define the 'spatial' energy flows between 
the lattice sites (n) and (n-|- 1) of the FF and SH compo- 
nents, respectively, and Jp(n) = — 27lm(y4^_B,*) deflnes 
the 'spectral' energy flow between FF and SH compo- 
nents due to quadratic nonlinearity. 

We seek vortex solutions where the energy flows in 
spatial and spectral domains form closed loops, such 
that there is no overall energy redistribution between 
the waveguides along the propagation direction. Such 
states correspond to stationary solutions for Eqs. ([T]) 
in the form An = Ao(n) exp [i6z -I- i(^i(n)], and Bn = 
BQ{n) ex.'p^ihz -\r iip2{n)\, where b is the real propaga- 
tion constant. After substituting these expressions in 
Eqs. ll]), wc obtained a set of stationary nonlinear equa- 
tions for the real amplitude [ylo(n) and i?o(n)] and 
phase [1^91.2 (n)] profiles of FF and SH components, which 
were solved with the numerical relaxation technique. The 
power flows for stationary solutions are expressed as: 
Ji{n,n+\) = 2ciAo(n)Ao(n+ 1) sin[i^i(n-fl) - ^i{n)\, 
J2{n,n+l) = 2c2i?o('^)-Bo(n- + l)sin[(^2(n- + l)-V52(?^)], 
and Jp{n) = 2'y\Ao{ri)\'^ BQ{n) sui[ip2{n) — 2(^i(n)]. Note 
that the power flows are non-zero only when there ap- 
pears nontrivial phase difference (not or tt) at adja- 
cent lattice sites for FF and SH fields. Whereas com- 
plex solutions were identified for quadratic couplers [7], 
only in-phase or twisted localized modes with vanish- 
ing power flows were previously found in quadratic lat- 
tices [3,4,6,8]. Wc also note that complex phase-twisted 
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states were found within the framework of single ex- 
tended discrete nonhnear Schrodinger equation [9] . 

We find new classes of solutions featuring nontrivial 
phase shift between the neighboring lattice sites. Some 
illustrative examples of nontrivial phase-twisted modes 
are presented in Fig. [Jl from which one can sec that the 
stationary solutions arc complex, in sharp contrast with 
the trivial phase twisted modes which have real trans- 
verse profiles. Thus we identify them as self-localized 
optical vortices, the spatial- spectral vortices (SSV). For 
the fundamental vortex solution presented in Fig. [TJa), 
the structure of energy flows can be qualitatively ex- 
plained as follows: (i) the energy flow of the FF wave 
is directed in one spatial direction (shown by top ar- 
row); (ii) the energy flow is oppositely directed at the SH 
wave (shown by bottom arrow). One can sec clearly that 
the energy flows between the fundamental and second- 
harmonic waves are counter-directed at different spatial 
locations (shown by vertical arrows). We also find a rich 
variety of different vortex states. The solution presented 
in Fig. HKb) has an asymmetric shape with one loop of en- 
ergy flow, which can be viewed as asymmetric SSV. The 
mode shown in Fig. [ijc) features two loops of the energy 
flow. We find that the SSV presented in Figs.[lja) and (c) 
have similar properties, whereas asymmetric SSVs are 
found to be unstable. Therefore, below we present the 
comprehensive analysis of existence and stability prop- 
erties of SSV with a single energy flow loop, as shown in 
Fig. Ha). 

For a given waveguide array characterized by specific 
values of coefficients in model Eqs. ([T|), the vortex so- 
lutions form families parametrized by the propagation 
constant b. An important characteristic of SSV is the to- 
tal power, as shown in Fig. [21(a). There exists a cutoff 
value of propagation constant b, which depends on the 
sign and absolute value of phase mismatch parameter, 
and the coupling coefficients. We also find that the power 
is a nonmonotonic function of the propagation constant, 
and there is a narrow region near cutoff (not even visible 
in the figure) where dl/db < 0. Our calculations show 
that the SSV have the power which is approximately two 
times higher than the power of fundamental (single-site) 
discrete solitons. Therefore such SSV can be considered 
as nontrivial bound states of fundamental solitons. 

The key point of this work is that SSV appear due to 
the simultaneous effects of parametric wave mixing and 
linear mode coupling between the waveguides. We plot 
the dependence of soliton power on the ratio of coupling 
coefficients for a fixed propagation constant and phase 
mismatch parameters in Fig. El^b). In this plot the solid 
and dotted lines corresponds to SSV with complex phase 
profiles, and we see that the existence of SSV requires 
that the coupling strength for the SH field exceeds a 
certain threshold. The SSV are found to bifurcate from 
the twisted soliton having or tt phase jump between 
lattice sites (shown with dashed line), above a critical 
value of the SH coupling coefficient. 

It is important to determine the stability proper- 
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Fig. 1. (Color online) Examples of SSV: column (a) sym- 
metric vortex; column (b) asymmetric single-charge vor- 
tex, and column (c) vortex with two power fiow loops. 
Here the first and second rows show the intensity and 
phase distributions of SSV, where lines with white and 
black circles show the FF and SH fields, respectively. In 
the third row, a scheme illustrates the energy flow of 
SSV. Here /3 = 0, 6 = 4, and C2 = 2ci = 1. 



ties of vortices. We perform the linear stability analy- 
sis by considering the perturbed solutions in the form 
Aniz) = {Ao(n)e'''^i(") -f /l(n)e*^}e*^^ and B„(z) = 
{Bo(n)e*'^2(") + f2{n)e^^^} e^'''\ Here S is the instabil- 
ity growth rate, and fi^2{n) are the corresponding mode 
profiles. These were determined by numerical solutions of 
the eigenvalue problem, obtained by linearizing Eqs. (P) 
with respect to small-amplitude perturbations. We indi- 
cate stable SSV with solid lines, and unstable SSV with 
dotted lines in Figs.[2l^a) and[2lb), and insets in these fig- 
ures show the real part of the perturbation growth rate. 
Our results show that SSV are stable when their power 
exceeds a certain threshold [see Fig. [21[a)], as the growth 
rate goes to zero. When the power is below a critical 
value, SSV becomes unstable and suffer from oscillatory 
instability (leading to its decay into the fundamental soli- 
tons). The variation of phase mismatch conditions does 
not change dramatically the stability-instability scenar- 
ios for SSV, and the whole instability domain shown in 
Fig.[2jd) confirms that SSV are stable when their power 
is above a certain value even for different phase mis- 
match conditions. The coupling strength also affects the 
stability of SSV. Fig. El^b) shows that for a fixed non- 
linear wave number shift the SSV become more stable 
for smaller coupling strength. Eventually when the cou- 
pling strength is below the critical value SSV transform 
into trivial phase-twisted discrete solitons (shown with 
dashed line), which are stable in the same parameter 
regime. Note that the instability domain of SSV increases 
with the increase of the SH coupling for the fixed cou- 
pling strength of FF field [Fig.[2fc)]. The results of linear 
stability analysis have been confirmed by the extensive 
numerical simulations of Eq. ([T]) using the beam propa- 
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Fig. 2. Existence and stability of SSV as shown in 
Fig. [BJa) : (a,b) Power and perturbation growth rate (in- 
set) vs. (a) the propagation constant at C2 = 1, /3 = 
and (b) the coupUng coefficient of SH field at 5 = 4, 
/? = 0. Solid lines correspond to stable SSV, dotted - un- 
stable SSV, and dashed - twisted solitons without vortex 
flows. (c,d) Stable (gray) and unstable (light gray) do- 
mains versus (c) the coupling coefficient of SH field at 
[3 = 0, and (d) the phase mismatch at C2 = f . In all the 
cases ci = 0.5. 

gation method. 

To verify that the SSV can be observed in experiment, 
we simulate their generation under practical conditions. 
We consider a 7cm long waveguide array formed by ti- 
tanium in-diffusion into the Z-cut LINBO3 surface, and 
the interchannel half-beat coupling length for FF wave is 
Lc = 22mm, which results in the coupling strength for FF 
wave of ci = 0.5. The phase-mismatch can be changed 
by tuning the sample temperature T, which is given by 
P = -AkL = 8.1(234 - TpC]). Additionally, strong 
linear coupling between the SH waveguide modes is re- 
quired (we assume C2 = 2ci). Although in previous ex- 
periments [10] C2 <C ci, strong coupling can be achieved 
for higher-order SH guided modes, which can be selec- 
tively excited when their propagation constants satisfy 
the phase-matching condition [11]. We consider FF wave 
as an input, and launch a Gaussian beams with a phase 
tilt along the sites, thus the input FF beam is in the form 
An{z = 0) = aexp[— ((x„ -I- a;o)/ro)^]exp(iQ;7T,), where a 
being amplitude, xq = 0.5 is the transverse dislocation, 
and a being the phase tilt. Fig. [3] shows an illustrative 
example for the excitation of SSV. Indeed, even as the 
stationary state is reached, the power flows between the 
lattice sites persist. The estimated input power required 
for the vortex generation is 2KW (corresponds to a = 3.5, 
ro = 3.0, and a = 7r/4). 

In conclusion, we have shown that spatial-spectral vor- 
tices can appear in quadratic nonlinear waveguide ar- 
rays. Such modes have a nontrivial phase between neigh- 
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Fig. 3. (Color online) Generation of spatial-spectral vor- 
tex with FF beam input. Top row: evolution of amplitude 
modulus for (a) FF and (b) SH components. Bottom row: 
evolution of energy flows for (c) FF and (d) SH compo- 
nents. The insets with circles show the output profiles. 
Here P = 0, and C2 = 2ci = 1. 

boring sites, facilitating non- vanishing energy flows along 
closed loops. We stress that the existence of SSV is only 
allowed due to the effect of parametric wave mixing. 
These results suggest new opportunities in using non- 
linear processes to control the flow of light. 

The authors acknowledge fruitful discussions with 
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the Australian Research Council. 

References 

1. F. Lederer, G. Stegeman, D.N. Christodoulides, G. As- 
santo, M. Segev, and Y. Silberberg, Phys. Rep. 463, 1 
(2008). 

2. T. Pertsch, U. Pescliel, and F. Lederer, Opt. Lett. 28, 
102 (2003). 

3. S. Darmanyan, A. Kobyakov, and F. Lederer, Phys. Rev. 
E 57, 2344 (1998). 

4. B. A. Malomed, P.G. Kevrekidis, D.J. Frantzeskakis, 
H.E. Nistazakis, and A.N. Yannacopoulos, Phys. Rev. 
E 65, 056606 (2002). 

5. Y. V. Kartashov, L. Torner, and V.A. Vysloukh, Opt. 
Lett. 29, 1117 (2004). 

6. T. Peschel, U. Peschel, and F. Lederer, Phys. Rev. E 57, 
1127 (1998). 

7. O. Bang, P.L. Christiansen, and C.B. Clausen, Phys. 
Rev. E 56, 7257 (1997). 

8. A. A. Sukhorukov, Yu.S. Kivshar, O. Bang, and CM. 
Soukoulis, Phys. Rev. E 63, 016615 (2001). 

9. M. Oster and M. Johansson, Phys. Rev. E 71, 025601 
(2005). 

10. R. Iwanow, R. Schiek, G.I. Stegeman, T. Pertsch, F. 
Lederer, Y. Min, and W. Sohler, Phys. Rev. Lett. 93, 
113902 (2004). 

11. A. Amoroso, A. Di Falco, G. Leo, G. Assanto, A. Parisi, 
A. Cine, and S.R. Sanseverino, IEEE Photonics Technol. 



3 



Lett. 15, 443 (2003). 



4 



